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We present results of a search for proton decays, p vK^ , using data from a 33 kton-year 
exposure of the Super-Kamiokande detector. Two decay modes of the kaon, — > fL^Vfi and K'^ 
tt^tt", were studied. The data were consistent with the background expected from atmospheric 
neutrinos; therefore a lower limit on the partial lifetime of the proton t /B{p vK^) was found to 
be 6.7 X 10^^ years at 90% confidence level. 

11.30. Fs, 11. 30. Pb, 13.30. Ce,14. 20. Dh,29.40.Ka 

One of the most unique predictions of Grand Unified Theories (GUTs) is baryon number violation. The minimal 
SU(5) GUT j|] predicts the dominant decay mode to be p — > e^Tr" with a predicted lifetime shorter than ^ 10^^ years. 
The current experimental lower limit is about 100 times longer than this Furthermore, the weak mixing angle 
predicted by this model does not agree with the experimental value and the three running coupling constants of the 
strong and electroweak forces do not meet exactly at a single point [||. Alternatively, the minimal supersymmctric 
(SUSY) SU(5) GUT makes a prediction for the weak mixing angle which is much closer to experimental results and 
predicts the lifetime of the proton decay into c^tt" to be more than four orders of magnitude longer than that in 
non-SUSY minimal SU(5) GUT Q. The minimal SUSY SU(5) model predicts the proton decay mode p —> vK^ to 
be dominant with the partial lifetime prediction varying from 0(10^^) to 0(10^^) yr M. 
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In this letter, we report the result of the search for proton decay through the channel p — > PK~^ in 535 live-days of 
data in Super-Kamiokande, corresponding to a 33 kt-year exposure. 

The Super-Kamiokande detector is a ring imaging water Cherenkov detector located in Kamioka Observatory, 
ICRR, Univ. of Tokyo, 1000 m (2700 m water equivalent) below the peak of Mt. Ikenoyama near Kamioka, Japan. 
50 kilotons of ultra pure water are held within a stainless steel tank of height 41.4 m and diameter 39.3 m. The 
tank is optically separated into two regions, the inner and outer detectors, by a photomultiplier tube (PMT) support 
structure and a pair of opaque plastic sheets. The inner detector has a height of 36.2 m and a diameter of 33.8 m. 
The outer detector completely surrounds the inner detector and is used to identify incoming and outgoing particles. 
On the wall of the inner detector, there are 11146 50 cm inward facing PMTs which cover 40% of the surface. The 
outer detector is lined with 1885 20 cm outward facing PMTs equipped with 60 cm x 60 cm wavelength-shifter plates 
to increase the collection efficiency of Cherenkov photons. 

In the search for rare proton decays, neutrinos produced by cosmic ray interactions in the upper atmosphere (at- 
mospheric neutrinos) typically represent the limiting background. Super-Kamiokande detects about 8 fully-contained 
atmospheric neutrino events per day These events have vertices which were reconstructed inside the 22.5 kiloton 
fiducial volume of the detector, defined to be 2 m from the PMT support structure. All the events were required to 
have no coincident light in the outer detector. 

After the event selection, events were reconstructed using PMT pulse height and timing information to determine 
the vertex position, number of Cherenkov rings, particle type, momentum, and number of decay electrons. A particle 
is classified as a showering(e-like) or a non-showering(/i-like) type Details of the detector, event selection, and 
event reconstruction are described in and 

The absolute scale of the momentum reconstruction was checked with several calibration sources such as the 
electrons from a LIN AC ^] , decay electrons from cosmic-ray muons which stopped in the detector volume (stopping 
muons), dE/dx of stopping muons, and the reconstructed mass of tt'^s produced by atmospheric neutrino interactions. 
The error in the absolute energy scale was estimated to be less than ±2.5%. The time variation of the energy scale 
was checked with muon decay electrons and was found to vary by ±1% over the exposure period. 

The momentum of the from p — > DK~^ is 340 MeV/c and is below the threshold momentum for producing 
Cherenkov light in water. Candidate events for this decay mode are therefore identified through the decay products 
of the K~^. Due to the smallncss of the hadronic cross-section of low- momentum K'^, we calculated that 97% of the 
K'^ exit from the ^®0 nucleus without interaction, and that 90% of these decay at rest. In this paper, we therefore 
searched for decays at rest through two dominant modes: K'^ M^^'m ^^^"^ ^ tt+tt''. Two separate methods 
were used to search for K'^ p^v^- 

We now describe the first method to search for p — > vK^ ; /^^^'^t • The momentum of the from the 

decay of the stopped is 236 MeV/c. If a proton in the P3/2 state of ^^O decays, the remaining ^^N nucleus is left 
in an excited state. This state quickly decays, emitting a prompt 6.3 MeV 7-ray. The signal from the decay particles 
of the should be delayed relative to that of the 7-ray due to the lifetime of the (r = 12 ns). The probability 
of a 6.3 MeV 7-ray emission in a proton decay in oxygen is estimated to be 41% [Q. By requiring this prompt 7-ray, 
almost all of the background events are eliminated. 

To determine the signal region and to estimate the detection efficiency, a total of 2000 p vK^ ; M^^'/x 
Monte Carlo (MC) events were generated. The same reduction and reconstruction as for the data were applied to 
these events. For these proton decay MC events, 96% were identified as having one ring. The resolution of the vertex 
fitting was 52 cm and the probability of misidentification of the particle type was 2.9%. Based on this simulation, 
a signal region of the reconstructed momentum of the /i+ was defined to be between 215 MeV/c and 260 MeV/c. 
Candidate events for this analysis were then selected by the following criteria: Al) one ring A2) /i-like A3) with one 
decay electron, A4) 215 MeV/c < momentum(/x) < 260 MeV/c, A5) Goodness of fit requirement, and A6) detection 
of a prompt 7-ray. The A5 criterion required a successful vertex fitting for the correct identification of the prompt 
7-ray. The prompt 7-ray signal was selected using the number of hit PMTs within a 12ns timing window which slides 
between 12ns and 120ns before the fj, signal. Criterion A6 required the number of hit PMTs to be more than 7. The 
total detection efficiency, including the selection criteria AI-A6, the branching ratio of K'^ l^^v^ (63.5%), the 
ratio of bound protons in oxygen to all of the protons in II2O (80%), and the emission probability of the 6.3 MeV 
7-ray (41%), was estimated to be 4.4%. The reason for the relatively low efhciency was due to the inefficiency of 
the detection of the prompt 7-rays. The number of background events was estimated using a 225 kt-year equivalent 
sample of atmospheric neutrino MC events. When estimating this background, the flux was normalized based on the 
observed deficit of f^'s. The normalization factors were 0.74 for charged current events and 1.17 for any other 
neutrino induced events The background contamination was estimated to be 0.4 events/33 kt-yr. When the same 
criteria were applied to the real data, no events were observed. Figure ^ shows the number of hit PMTs within the 
timing region of the prompt 7-ray for proton decay MC events, simulated atmospheric neutrino events, and real data. 
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The partial lifetime (t/B) was calculated: 



t/(3 > (A X eB„,)/N^and, 



(1) 



where A is the exposure in proton-years, cBm is the detection efficiency (4.4%) and Ncand is the 90% Poisson upper 
limit on the number of candidate events (2.3). The lower limit of the partial lifetime obtained for this mode was found 
to be 2.1 X 10^2 yr at 90% C.L. 

The overall detection efficiency of the prompt 7-ray tagging method is rather small. Therefore, as a second method 
we searched for an excess of events of mono-energetic 236 MeV/c /Lt+'s produced by the stopped decay for events 
with no observed prompt 7-ray. For this analysis, the selection criteria were: Al to A4 described above and A7 no 
prompt gamma-ray signal. The detection efficiency of this mode (including the branching ratio of —>■ i-i^i'fi ) was 
estimated to be 40%. 

To estimate the excess of proton-decay signal, the number of events in three momentum regions, 200 MeV/c to 215 
MeV/c, 215 MeV/c to 260 MeV/c, and 260 MeV/c to 300 MeV/c, were summed separately for p vK+ \K+ -> ^jl+v^ 
MC, atmospheric v MC, and data. The method was applied to fit parameters. The function was defined as 
follows: 



where a and b are the fit parameters, N f N f ,Nf^'' are the numbers of events of real data, proton decay MC, 
and atmospheric v MC, respectively, in each momentum region i. 

The minimum x^, Xmin = O-^i was in the unphysical region [b ■ TVl"''^^ — —16.3). The Xmin ™ ^^le physical region 
(6 • N^'^'^^ = 0), XminiP^y^) 2.0. The data were therefore consistent with no excess of p — ^ events. The 

90% C.L. upper limit on the number of proton decay events(iVcan!i) was obtained by requiring x^ — XminiP^V^) — 3-7 
based on the prescription described in Ref. H] . The 90% upper limit on the number of candidates was estimated to 
be 13.3. The momentum distribution of the events which satisfy criteria Al to A3 is shown in Figure |[ Also shown is 
the expected signal of proton decay at the 90% C.L. The lower limit of the partial lifetime oi p DK^; l^^^tj, 
using the above method was found to be 3.3 x lO'^^ years at 90% C.L. 

Finally, we describe the search for p — > vK^ ; ^ tt+tt" . The tt" and 7r+ from the decay at rest have 
equal and opposite momenta of approximately 205 MeV/c. The decay 7's from the tt" reconstruct this momentum. 
The TT"*", barely over Cherenkov threshold with (3 ~ .86, emits very little Cherenkov radiation. However it does decay 
into a muon (tt^ M^'^m) which decays into a positron (/i^ — > e^VeVfj). Detection of this positron is possible. 
Furthermore, a small amount of Cherenkov radiation from the 7r+, sometimes visible as a collapsed Cherenkov ring, 
can be detected in the direction opposite that of the 7r°. To quantify this, we defined "backwards charge" {Qb) as the 
sum of the photoelectrons detected by the PMTs which lay within a 40° cone whose axis was the opposite direction of 
the reconstructed direction of the tt". This charge was corrected for light attenuation in the water, angular dependence 
of photon acceptance, and photocathode coverage. 

To determine the signal region and estimate the detection efficiency, a total of 1000 p vK^ ; tt+tt*' MC 

events were generated, 771 of which had vertices which reconstructed inside the fiducial volume of the detector. For 
these 205 MeV/c 7r° events, the resolution of vertex fitting was 29 cm and 66% of the events were identified as 2-ring 
events. 

The selection criteria for this type of event were defined: Bl) 2 e-like rings, B2) with 1 decay electron, B3) 85 
MeV/c^ < mass^.., < 185 MeV/c^, B4) 175 MeV/c^ < momentum^^ < 250 MeV/c^ and B5) 40 p.e. < Qb < 100 
p.e. Criteria B1,B3, and B4 required the with the monochromatic momentum expected. Criterion B2 required 
the decay of the 7r+ into muon into positron. Criterion B5 required the Cherenkov light from the 7r+. The tt" mass 
resolution was determined to be 135±21 MeV/c^. By passing the proton decay MC events through these selection 
criteria, the detection efficiency was determined to be 31%. The largest contribution to the inefficiency was the 
inefficiency of detection of two 7-rays from the decay of the tt*'. Including the kaon branching ratio of 21.2% into 
tt+tt", the total detection efficiency for this mode was estimated to be 6.5%. 

Charged current interactions such as ly^^N —> /iiW^ from atmospheric neutrinos can imitate a kaon decay mode of 
this type. The selection criteria were applied to the sample of atmospheric neutrino MC. The flux was normalized 
in the same manner as in the prompt 7-ray search. The number of background events expected was estimated to be 
0.7 events/33kt-yr. Figures ^ and b show IpVol vs. backwards charge (Qf,) for proton decay MC and atmospheric 
neutrino MC, respectively. 

The selection criteria for this decay mode were applied to the data. Figure ^ shows the results of the final two 
cuts. No events passed. The single event which lay close to the cuts was examined visually with an event display. 
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We found no additional evidence that it could be a signal event that fell outside of the cuts. The backwards charge 
appeared to be a fragment of one of the rings and not from a small collapsed ring. Based on these numbers, the lower 
limit of the partial lifetime of tt+tt" was estimated to be 3.1 x 10^^ yr at the 90% C.L. based on Eqn. |^. 

The combined 90% C.L. upper limit for the number of proton decay candidates (xumit) was calculated by integrating 
the likelihood function to the 90% probability level: 



0.90, (3) 



iV.(x)^Xf'^+ X, (4) 

where P{N, x) is the probability function of Poisson statistics, N°^^ is the number of observed candidates, X^"^ is 
the number of estimated background events, ei is the detection efficiency, and -B^ is the meson branching ratio. The 
index i stands for the z-th method. The lower limit of partial lifetime was calculated with 

1 ^ 

r/B^ {J2e,B''JA, (5) 

Xlimit . T 

where A is the exposure in proton-years. The combined lower limit of the partial lifetime for t/B(jj DK^) using 
the three independent methods was 7.3 x 10'^^ yr at the 90% C.L. 

The main sources of systematic errors in the calculation of the lifetime limits were: 1) uncertainty in the energy 
calibration, 2) uncertainty in the detection efficiency of the decay electron from the stopped muon, 3) uncertainty in 
the atmospheric neutrino fluxes and interaction cross sections used in the Monte Carlo simulation, 4) uncertainty in 
the emission probability of prompt 7-ray from ^^O. 

For the K+ — > TT^TT mode search, the detection efficiency changed by less than ±1% if the reconstructed momentum 
criterion was shifted by ±2.5%. For the —* ji^v^ mode, if the reconstructed momentum was shifted by ±2.5%, 
the lifetime limit obtained from this analysis was changed by ±10% at most. The systematic error of the detection 
efficiency of decay electrons from muons was estimated to be 1.5% by comparing the fraction of cosmic-ray muon 
events with decay electrons between MC simulation and real data. The uncertainty of the estimation of number of 
background events came from the uncertainty of the neutrino flux and the interaction cross sections. This uncertainty 
was estimated by comparing the number of events of the normalized atmospheric neutrino background with the real 
data at each reduction step. For the prompt 7-ray tagging method and K'^ — > tt+tt" mode search, the data and 
atmospheric v Monte Carlo agreed within statistical errors at each step. In searching for the excess of 236 MeV/c 
/Lt+s, the shape of the momentum distribution and the absolute normalization was left as a free parameter in the 
fitting. Therefore, the uncertainties of neutrino flux and the interaction cross-sections did not significantly affect the 
lower limit of the partial lifetime. The systematic error in the 6.3 MeV 7-ray emission probability was estimated to 
be ±15% 1^. This error directly affected the detection efficiency of the prompt 7-ray tagging method. Also it is 
expected that there are some other promt 7-rays, whose energies are 9.93 MeV, 7.03 MeV, 7.01 MeV 0. Including 
these deexcitation modes, the detection efficiency of the prompt 7-ray tagging method should be higher. However, 
these emission probabilities are more than 10 times smaller than the probability of that of 6.3 MeV 7-ray and they 
have much larger uncertainties in the emission probability, which are estimated to be larger than 30% [|j . Therefore, 
we neglected these modes in the analysis. However, these modes were taken into account in the systematic uncertainty 
of the lifetime limit, which was estimated to be ^3%. 

Considering all of these effects, the systematic error of the lower limit of the partial lifetime was estimated to be 
lg°%. As a result, the lower limit of the partial lifetime for p — > vK^ was estimated to be 6.7 x 10'^^yr(90%C.L.). 

In this paper, we have reported the result of a search for proton decay into vK^ in a 33 kt-year exposure of the 
Super-Kamiokande detector. The data are consistent with the background expected from atmospheric neutrinos and 
no evidence for proton decay was observed. We set the lower limit of partial lifetime for p vK'^ to be 6.7 x 10^^ 
yr at 90% C.L. This limit is more than six times longer than the previously published best limit (1.0 xlO'^^ yr.(90% 
C.L.)) and will help to constrain SUSY GUT models. 

We gratefully acknowledge the cooperation of the Kamioka Mining and Smelting Company. The Super-Kamiokande 
experiment was built from and has been operated with, funding by the Japanese Ministry of Education, Science, Sports 
and Culture, and the United States Department of Energy. 
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FIG. 1. Number of hit PMTs within the timing region of the prompt 7-ray search in p —+ i/K ,K ^ Pi_i. Histogram 
shows the proton decay Monte Carlo with 6.3MeV 7-ray, shaded histogram shows 33kt-yr equivalent atmospheric neutrino 
Monte Carlo events, and data points with error bars show 33kt-yr data of Super-Kamiokande. For the proton decay Monte 
Carlo, t/B(p vK*) of 2.1 x 10^^ yr is assumed. 
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FIG. 2. Reconstructed momentum distribution. Solid (dotted) line shows the estimated 90% C.L. number of proton decay 
{t/B{p — » vK'^) = 3.3 X 10^^ yr) + atmospheric v (atmospheric v) events; the black points with error bars show the data with 
the statistical errors. 
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(b) Atmospheric v MC 
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(c) 33 kt-yr Super-K data 
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FIG. 3. Backwards charge versus iv^ momentum for (a) p — > i'K'^;K^ 
atmospheric u Monte Carlo, and (c) 33 kt-yr data from Super-Kamiokande. 
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